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Abstract 

Localized patterns and nonlinear oscillation formation on the bounded free 
surface of an ideal incompressible liquid are analytically investigated . Cnoidal 
modes, solitons and compactons, as traveling non-axially symmetric shapes are 
discused. A finite-difference differential generalized Korteweg-de Vries equation is 
shown to describe the three-dimensional motion of the fluid surface and the limit 
of long and shallow channels one reobtains the well known KdV equation. A tenta- 
tive expansion formula for the representation of the general solution of a nonlinear 
equation, for given initial condition is introduced on a graphical-algebraic basis. 
The model is useful in multilayer fluid dynamics, cluster formation, and nuclear 
physics since, up to an overall scale, these systems display liquid free surface 
behavior. 

PACS : 47.55.Dz, 68.10.Cr, 47.20.Ky, 47.20.Dr, 97.60.j, 83.70, 36.40.-c 
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1 Introduction 



Liquid oscillations on bounded surfaces have been studied intensively, both theoretically 
[1-3] and experimentally [4-6]. The small- amplitude oscillations of incompressible drops 
maintained by surface tension are usually characterized by their fundamental linear 
modes of motion in terms of spherical harmonics [1-3]. Nonlinear oscillations of a liquid 
drop introduce new phenomena and more complicated patterns (higher resonances, 
solitons, compactons, breakup and fragmentation, fractal structures, superdeformed 
shapes) than can be described by a linear theory. Nonlinearities in the description 
of an ideal drop demonstrating irrotational flow arise from Bernoulli's equation for 
the pressure field and from the kinematic surface boundary conditions [7]. Computer 
simulations have been carried for non-linear axial oscillations and they are in very good 
agreement with experiments [4-6]. 

The majority of experiments show a rich variety of complicated shapes, many re- 
lated to the spinning, breaking, fission and fusion of liquid drops. There are experiments 
[6] and numerical simulations [2] where special rotational patterns of circulation emerge: 
a running wave originates on the surface of the drop and then propagates inward. Re- 
cent results (superconductors [8], catalytic patterns [9], quasi-molecular spectra [10], 
numerical tests on higher order non-linear equations [11] and analytical calculations 
on the non-compact real axis [12-13]) show shape-stable traveling waves for nonlinear 
systems with compact geometry. Recent studies showed that a similar one-dimensional 
analysis for the process of cluster emission from heavy nuclei and quasi-molecular spec- 
tra of nuclear molecules yields good agreement with experiment [10]. Such solutions are 
stable and express to a good extent the formation and stability of patterns, clusters, 
droplets, etc. However, even localised, they have nor compact support neither period- 
icity (excepting some intermediate steps of the cnoidal solutions, [10,13]), creating thus 
difficulties when analysing on compact surfaces. 

In the present paper we comment on the cnoidal-towards-solitons solution inves- 
tigated in [9,12], especially from the energy point of view. We introduce here a new 
nonlinear 3-dimensional dynamical model of the surface, in compact geometry (pools, 
droplets, bubbles, shells), inspired by [12], and we investigate the possibilities to obtain 
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compacton-like solutions for this model. We also study the scale symmetries of such 
solutions. 

The model in [9,12] consider the nolinear hydrodynamic equations of the surface of 
a liquid drop and show their direct connection to KdV or MKdV systems. Traveling so- 
lutions that are cnoidal waves are obtained [10,13] and they generate multiscale patterns 
ranging from small harmonic oscillations (linearized model), to nonlinear oscillations, 
up to solitary waves. These non-axis-symmetric localized shapes are described by a 
KdV Hamiltonian system, too, which results as the second order approximation of the 
general Hamiltonian, next corrextion from the linear harmonic shape oscillations. Such 
rotons were observed experimentally when the shape oscillations of a droplet became 
nonlinear [4,6,8,13]. 

2 Liquid drop cnoidal and soliton solutions from 
Hamiltonian approach 

The dynamics governing one-dimensional surface oscillations of a perfect (p =const.), 
irrotational fluid drop (or bubble, shell) can be described by the velocity field <J> and a 
corresponding Hamiltonian [1-3,7,10,13]. By expanding the Hamiltonian and dynamical 
equations in terms of a small parameter, i.e. the amplitude of the perturbation r\ over 
the radius of drop R , the usual linear theory is recovered in the first order. Higher 
order non-linear terms introduce deviations and produce large surface oscillations like 
cnoidal waves [7]. These oscillations, under conditions of a rigid core of radius Rq — h 
and non-zero angular momentum, transform into solitary waves. In the following, by 
using the calculation developed in [10], we present the Hamiltonian approach for the 
liquid drops nonlinear oscillations. However, this approach is different from the nuclear 
liquid drop model point of view in [10], since we do not use here the nuclear interaction 
(shell corrections) responsible for the formation of different potential valleys. 

The total hydrodynamic energy E consists of the sum of the kinetic T and po- 
tential U energies of the liquid drop. The shape function is assumed to factorize, 
r(0,(f),t) = Rq(1 + g(9)r)(<f>,t)). All terms that depend on 6 are absorbed in the coeffi- 
cients of some integrals and the energy reduces to a functional of i] only. The potential 
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energy is given by the surface energy Us — cr(A v — Ao)\v Q , where a is the surface pres- 
sure coefficient, A v is the area of the deformed drop, and Ao the area of the spherical 
drop, of constant volume V . The kinetic energy T = p §^ $V$ • dS/2, [1-3,10,13], the 
kinematic free surface boundary condition $ r = d t r + (d e r)& e /r 2 + {d^r^^/r 2 sin 
and the boundary condition for the radial velocity on the inner surface d r §\ r= R -h = 0, 
[7], result in the expression [2,3,10] 

T= R%P r« j*« Ro^Vt sin9 + ^grj^jl - sing) 

2 Jo Jo sji + gW + 9% 

If the total energy, written in the second order in 77, is taken to be a Hamiltonian 
H[i]], the time derivative of any quantity ^[77] is given by F t = [F,H]. Defining F = 
J 2?r r]((j) — Vt)d(f) it results ([10], last reference) 



dF r 27T r 2lT 

— = / i] t d(j) = / (2C 2 m + 6C 3 i]m - 2C 4 ?7000)(i0 = 0, (2) 
at Jo Jo 



which leads to the KdV equation. Here C 2 = <rH^{Sl$ + Sq'°/2) + B%pV 2 C%'~l/2, 
C 3 = aR 2 Sl'° 2 /2 + R 6 pV 2 (2S 3 !l 1 2 Ro + S 5 J 2 2 3 + R S 6 J 2 2 3 ) /2, C 4 = <jR 2 S^'°/2, with S*j = 
Rq 1 Jq h l g t glsm k 6d6. Terms proportional to r\rf^ can be neglected since they introduce a 
factor rjl/L 2 which is small compared to 77$, i.e. it is in the third order. In order to verify 
the correctness of the above approximations, we present, for a typical soliton solution 
r)(<j),t), some terms occuring in the expresion of E, Fig. 1. All details of calculation 
are given in [10,13]. Therefore, the energy of the non-linear liquid drop model can be 
interpreted as the Hamiltonian of the one-dimensional KdV equation. The coefficients 
in eq.(2) depend on two stationary functions of 6 (the depth h{6) and the transversal 
profile g(0)), hence, under the integration, they involve only a parametric dependence. 

The KdV equation has the following cnoidal wave (Jacobi elliptic function) as 
exact solution 



77 = a 3 + (a 2 - a 3 )sn 2 (^j 3 ^ ° 2 (0 - Vt) rnj , (3) 

where ai,a 2 ,a 3 are constants of integration, m 2 = (a 3 — a 2 )/(a 3 — cci). This solution 
oscillates between a 2 and 0:3, with a period T = 2K{m) \J ^ 3^ C3 ", where K(m) is the 
period of a Jacobi elliptic function sn(x\m). The parameter V is the velocity of the 
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cnoidal waves and ol\ + at 2 + «3 = 
approaches 



. In the limit a\ = a 2 = the solution eq.(3) 



r] = 7] sech 




(<j>-Vt) , 



(4) 



which is the soliton solution of amplitude r] . Small oscillation occur when a 3 — * a 2 and 
m — > 0,T — > 7r/2. Consequently, the system has two limiting solutions, a periodic and 
a localized traveling profile, which deform one into the other, by the initial conditions 
and the velocity parameter V . A figure showing the deformation from the / = 5 cnoidal 
mode towards a soliton is shown in Figs. 2. 

The cnoidal solution eq.(3) depends on the parameters ctj subjected to the volume 
conservation and the periodicity condition of the solution (for the final soliton state 
this condition should be taken as a quasi-periodicity realised by the rapidly decreasing 
profile. This a problem of the basic model, [10]). The periodicity restriction reads 



Hence, a single free parameter remains, which can be taken either one out of the three 
a's, V or i] . Equatorial cross-sections of the drop are shown in Fig. 2b for the cnoidal 
solution at several values of the parameter 770. All explicite calculations are presented 
in detail in [10]. 

In Fig. 3 we present the total energy plotted versus the parameters a±,a 2 for 
constant volume. From the small oscillation limit (a 2 — 3 in the figure) towards the 
solitary wave limit (a 2 = 1 in the figure) the energy increases and has a valley for 
ai ~ 0.1 and a 2 G (1.2, 1.75) (close to the I = 2 mode). In order to introduce more 
realistic results, the total hydrodynamic energy is plotted versus a±,a 2 for constant 
volume, too but we marked those special solutions fulfilling the periodicity condition. 
In Fig. 4 we present the total energy valley, from the small oscillations limit towards 
the solitary wave limit. We notice that the energy constantly increases but around 
a 2 G (1.2, 1.75) (close to the linear I = 2 mode) it has a valley providing some stability 
for solitary solution (also called roton [13]). 




n — 1,2, ... , 2y/a 3 — a[. 



(5) 
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3 The three-dimensional nonlinear model 



In the following we introduce a sort of generalized KdV equation for fluids. We consider 
the three-dimensional irrotational flow of an ideal incompressible fluid layer in a semi- 
finte rectangular channel subjected to uniform vertical gravitation (g in z direction) 
and to surface pressure [12]. The depth of the layer, when the fluid is at rest is z = 
h. Boundary conditions at the finite spaced walls consist in annilation of the normal 
velocity component, i.e. on the bottom of the layer (z = 0) and on the walls x = x ±L/2 
of the channel of width L. The following results remain valid if the walls expand 
arbitrary, e.g. L — > oo, and the flow is free. We choose for the potential of the 
velocities the form 

M kn(x - x ) V2kn(y - y ) knz 

9 = a k (t) cos cosh cos——, (6) 

k>o L L L 

where are arbitrary functions of time and L is a free parameter. Eq.(7) fulfils 

A$ = and the above boundary conditions at the walls. However there is another 
boundary condition at the free surface of the fluid [7] 

($ 2 -T] t - Vx^x)z=h+r, = 0, (7) 

where r/(x, t) describes the shape of the free surface. By introducing the function 

ft +\ sr- a k(t)kn{. k7T(x-x ) V2kn(y-y )\ f . 
f(x,t) = 2^ — - — (^sm cosh J, (8) 

the velocity field on the free surface can be written 

®x\z=h+ri = -cosh(zd x )f, 

-&z\z=h+r, = ~ smh(zd x )f. (9) 

Eqs (10) do not depend on L and the case L — > oo of unbounded channels and free 
travelling profiles remains equaly valid. Since the unique force field in the problem 
is potential, the dynamics is described by the Bernoulli equation, which, at the free 
surface, reads 

<$>Xt + + + 9Vx + ~Px = 0. (10) 
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Here P is the surface pressure obtained by equating P's first variation with the local 
mean curvature of the surface, under the restriction of the volume conservation 

OJ} xx 



z = h+ v (l + V^f 2 ' 

and a is the surface pressure coefficient. The pressure in eq.(12) approaches —<Ji] xx , for 
small enough relative amplitude of the deformation rj/h. In order to solve the system 
of the two partial differential equations (8,11) with respect to the unknown functions 
f(x,t) and r](x,t), we consider the approximation of small perturbations of the surface 
compared to the depth, a = max\r]^ k \x,t)\ « h, where k = 0, ...,3 are orders of 
differentiation. Inspired by [12] and using a sort of perturbation technique in a/h, we 
obtain from eqs.(6-ll) the generalised KdV equation 

Vt + ~r sm(hd)i] + -r{i]x cosh(hd)rj + rj cosh(hd)i] x )0. (12) 

I i I i 

If we approximate sin(/i<9) ~ hd — |(/i<9) 3 , cosh(/i<9) ~ 1 — ^(hd) 2 , we obtain, from 
eq.(9), the polynomial differential equation: 

_ o, ~~ ,~ h 3 ~ c 6 2 h 3 /„ _ _ \ n /io , 
ar} t + 2c e hr}r] x + c Q ehr] x - c e—r] xxx — \n x r] xx + i]r] xxx I =0, (13) 

where e = f • The first four terms in eq.(20) correspond to the zero order approximation 
terms in /, obtained from the boundary condition at the free surface, i.e. the traditional 
way of obtaining the KdV equation in shallow channels. 

In order to find an exact solution for eq.(12) we can write it in the form: 

Ahux(X) + + + ux{x) u(X + k) + u(X-k) 

2i 2 

2 

where X = x + Ac t and A is an arbitrary real constant. We want to stress here that 
eq.(14) is a finite-difference differential equation, which is rather the exception than 
the rule fir such systems. Hence, it may contain among its symmetries, the scaling 
symmetry. Actualy, the first derivative of u(X) is shown to be alinear combination of 
translated versions of the original function. In this way, the theory of such equations 



can be related with the wavelet, or other self-similarity systems, theory, [13]. In the 
following we study the solutions with a rapid decreasing at infinity and make a change 
of variable: v = e BX for x G (—00, 0) and v = e~ BX for x G (0, 00), with B an arbitrary 
constant. Writing u(X) = —hA + f(v), and choosing for the solution the form of a 
power series in v: 

00 

f(v) = J2 a nV n , (15) 
n=0 

we obtain a nonlinear recurrsion relation for the coefficients a n : 

sin(Bhk) \ 



\Ahk+ — \-—L\ ak 



= — ^ ^( cosh (Bh(k — n)) + cosh(Bh(k — 1)) \a n a k ^ n . (16) 

n=l ^ ' 

With the coefficients given in eq.(16) the general solution rj can be written analyticaly. 
In order to verify the consistency of this solution we study a limiting case of the relation, 
by replacing sin and cosh expressions with their lowest nonvanishing terms in their power 
expansions Thus, eq.(16) reduces to 

k-l 
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n ' = — 77 X! ""><">■• ( 17 ) 

L ) n=l 



B 2 h S k(k 2 



and 



fe-i 



= k < 18 » 

is the solution of the above recurrence relation. In this approximation, the solution of 
eq.(12) reads 

y ' £^ V J 2 (cosh( J BX/2)) 2 V ; 

which is just the single-soliton solution of the KdV equation and it was indeed obtained 
by assuming h small in the recurrence relation (16). Hence, we have shown that the KdV 
equation describing the shallow liquids can be generalised for any depths and lengths. 
This result may be the starting point to search for more interesting symmetries. It 
would be interesting to interpret the generalized-KdV eq.(12) as the Casimir element 
of a certain algebra. 



4 Compacton and self- similar solutions 



Eq.(12) has a special character, namely contains both infinitesimal and finite difference 
operators. This particularity relates it to another field of nonlinear systems, that is 
scaling functions and wavelet basis, functions or distributions with compact support and 
self-similarity properties. In the following we investigate a particular case of eq.(12), 
that is when h <C 77, h <C 5, where 5 is the half-width of the solution, if this has bounded 
or compact support. In this approximations, from eq.(12) we keep only the terms 

—in + Vx + \mx - ttVxVxx + \mx - TiWxxx + e> 3 - o. (20) 

c h 2 h 2 

This equation is related to another intergable system, namely the K(2,2) equation, 
investigated in [11] 

Vt + (V 2 )x + (V 2 )xxx = 0. (21) 

The main property of the K(2,2) equation is the equal occurence of non-linearity, dis- 
persion and the existence of a Lagrangian and Hamiltonian system associated with it. 
The special solutions of this equation are the compactons 

r/ c = ^cos 2 f^V \x- m t\>2„, (22) 
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and i] c = otherwise. This special solutions have compact support and special proper- 
ties concerning the scattering between different such solutions. As the authors comment 
in [11], the robustness of these solutions makes it clear that a new mechanism is un- 
derlying this system. In this respect, we would like to add that, taking into account 
eq.(12), this new mechanism might be related to selfsimilarity and multiscale properties 
of nonlinear systems. 



5 Conclusions 



In the present paper we introduced a non-linear hydrodynamic model describing new 
modes of motion of the free surface of a liquid. The total energy of this nonlinear liquid 
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drop model, subject to non-linear boundary conditions at the free surface and the inner 
surface of the fluid layer, gives the Hamiltonian of the Korteweg de Vries equation. We 
have studied the stability of the cnoidal wave and solitary wave solutions, from the 
point of view of minima of this Hamiltonian. 

The non-linear terms yield rotating steady-state solutions that are cnoidal waves 
on the surface of the drop, covering continuously the range from small harmonic oscilla- 
tions, to anharmonic oscillations, and up to solitary waves. The initial one-dimensional 
model [10] was extend to a three-dimensional model. A kind of new generalized KdV 
equation, together with some of its analytical solutions have been presented. We also 
found a connection between the obtained generalized KdV equation, and another one 
(i.e. K(2,2)), in a certain approximation. In this case, compacton solutions have been 
found and new symmetries (e.g. self-similarity) were put into evidence. 

The analytic solutions of the non-linear model presented in this paper, make pos- 
sible the study of clusterization as well as to explain or predict the existence of new 
strongly deformed shapes, or new patterns having compact support or finite wavelength. 
The model applies not only in fluid and rheology theories, but may provide insight into 
similar processes occurring in other fields and at other scales, such as the behavior 
of superdeformed nuclei, supernova, preformation of cluster in hydrodynamic models 
(metallic, molecular, nuclear), the fission of liquid drops (nuclear physics), inertial fu- 
sion, etc. 
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FIGURE CAPTIONS 



Fig. 1 

The order of smallness of four typical terms depending on <fi and occuring in the 
Hamiltonian, eqs.(l,2). Order zero holds for i] 2 , order 1 for 77^, order 2 for rf, and order 
3 for 7777^. 

Figs. 2 
2a. 

The transition of the cnoidal solution, from a I — 5 mode to the soliton limit: 
shape of the cross-section for 9 = n/2 function as a function of a 2 with fixed by 
the volume conservation and periodicity conditions. 

2b. 

Cnoidal solutions (cross-sections of El for 9 = 7r/2) subject to the volume con- 
servation constraint. Results for the / = 6 mode to the I = 2 mode and a soliton are 
shown. The corresponding linear modes, i.e. spherical harmonics, are superimposed on 
the non-linear solutions. 

2c. 

Pictorial view of a soliton deformation of a drop, on the top of the original unde- 
formed sphere. The supporting sphere for the soliton has smaller radius because of the 
volume conservation. 



Fig. 3 

The energy plotted versus a±,a 2 for constant volume. From the small oscillation 
limit (a 2 — 3) towards the solitary wave limit (a 2 = 1) the energy increases and has a 
valley for a,\ ~ 0.1 and a 2 G (1.2, 1.75) (close to the / = 2 mode). 
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Fig. 4 

The total energy plotted versus a±, a 2 for constant volume (small circles). Larger 
circles indicate the patterns fulfilling the periodicity condition. From the small oscilla- 
tions limit (a 2 — 3) towards the solitary wave limit (a 2 = 1) the energy increases but 
for a 2 G (1.2, 1.75) (close to I = 2 mode) it has a valley. 
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